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Abstract — We consider a general multiple antenna network 
with multiple sources, multiple destinations and multiple relays in 
terms of the diversity-multiplexing tradeoff (DMT). We examine 
several subcases of this most general problem taking into account 
the processing capability of the relays (half-duplex or full-duplex), 
and the network geometry (clustered or non-clustered). We first 
study the multiple antenna relay channel with a full-duplex relay 
to understand the effect of increased degrees of freedom in 
the direct link. We find DMT upper bounds and investigate 
the achievable performance of decode-and-forward (DF), and 
compress-and-forward (CF) protocols. Our results suggest that 
while DF is DMT optimal when all terminals have one antenna 
each, it may not maintain its good performance when the degrees 
of freedom in the direct link is increased, whereas CF continues 
to perform optimally. We also study the multiple antenna relay 
channel with a half-duplex relay. We show that the half-duplex 
DMT behavior can significantly be different from the full-duplex 
case. We find that CF is DMT optimal for half-duplex relaying 
as well, and is the first protocol known to achieve the half- 
duplex relay DMT. We next study the multiple-access relay 
channel (MARC) DMT. Finally, we investigate a system with 
a single source-destination pair and multiple relays, each node 
with a single antenna, and show that even under the idealistic 
assumption of full-duplex relays and a clustered network, this 
virtual multi-input multi-output (MIMO) system can never fully 
mimic a real MIMO DMT. For cooperative systems with multiple 
sources and multiple destinations the same limitation remains to 
be in effect. 

Index Terms — cooperation, diversity-multiplexing tradeoff, 
fading channels, multiple-input multiple-output (MIMO), relay 
channel, wireless networks. 



I. Introduction 

Next-generation wireless communication systems demand 
both high transmission rates and a quality-of-service guar- 
antee. This demand directly conflicts with the properties of 
the wireless medium. As a result of the scatterers in the 
environment and mobile terminals, signal components received 
over different propagation paths may add destructively or 
constructively and cause random fluctuations in the received 
signal strength [1]. This phenomena, which is called fading, 
degrades the system performance. Multi-input multi-output 
(MIMO) systems introduce spatial diversity to combat fading. 
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Additionally, taking advantage of the rich scattering environ- 
ment, MIMO increases spatial multiplexing [2], [3]. 

User cooperation/relaying is a practical alternative to MIMO 
when the size of the wireless device is limited. Similar to 
MIMO, cooperation among different users can increase the 
achievable rates and decrease susceptibility to channel varia- 
tions [4], [5]. In [6], the authors proposed relaying strategies 
that increase the system reliability. Although the capacity of 
the general relay channel problem has been unsolved for over 
thirty years [7], [8], the papers [4], [5] and [6] triggered a vast 
literature on cooperative wireless systems. Various relaying 
strategies and space-time code designs that increase diversity 
gains or achievable rates are studied in [9] -[36]. 

As opposed to the either/or approach of higher reliability or 
higher rate, the seminal paper [37] establishes the fundamental 
tradeoff between these two measures, reliability and rate, 
also known as the diversity-multiplexing tradeoff (DMT), for 
MIMO systems. At high SNR, the measure of reliability 
is the diversity gain, which shows how fast the probability 
of error decreases with increasing SNR. The multiplexing 
gain, on the other hand, describes how fast the actual rate 
of the system increases with SNR. DMT is a powerful tool 
to evaluate the performance of different multiple antenna 
schemes at high SNR; it is also a useful performance measure 
for cooperative/relay systems. On one hand it is easy enough 
to tackle, and on the other hand it is strong enough to show 
insightful comparisons among different relaying schemes. 
While the capacity of the relay channel is not known in 
general, it is possible to find relaying schemes that exhibit 
optimal DMT performance. Therefore, in this work we study 
cooperative/relaying systems from a DMT perspective. 

In a general cooperative/relaying network with multiple 
antenna nodes, some of the nodes are sources, some are 
destinations, and some are mere relays. Finding a complete 
DMT characterization of the most general network seems 
elusive at this time, we will highlight some of the challenges 
in the paper Therefore, we examine the following important 
subproblems of the most general network. 

• Problem 1: A single source-destination system, with one 
relay, each node has multiple antennas, 

• Problem 2: The multiple-access relay channel with mul- 
tiple sources, one destination and one relay, each node 
has multiple antennas, 

• Problem 3: A single source-destination system with mul- 
tiple relays, each node has a single antenna, 

• Problem 4: A multiple source-multiple destination sys- 
tem, each node has a single antenna. 

An important constraint is the processing capability of 
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the relay(s). We investigate cooperative/relaying systems and 
strategies under the full-duplex assumption, i.e. when wireless 
devices transmit and receive simultaneously, to highlight some 
of the fundamental properties and hmitations. Half-duplex sys- 
tems, where wireless devices cannot transmit and receive at the 
same time, are also of interest, as the half-duplex assumption 
more accurately models a practical system. Therefore, we 
study both full-duplex and half-duplex relays in the above 
network configurations. 

The channel model and relative node locations have an 
important effect on the DMT results that we provide in this 
paper In [38], we investigated Problem 3 from the diversity 
perspective only. We showed that in order to have maximal 
MIMO diversity gain, the relays should be clustered around 
the source and the destination evenly. In other words, half 
of the relays should be in close proximity to the source and 
the rest close to the destination so that they have a strong 
inter-user channel approximated as an additive white Gaussian 
noise (AWGN) channel. Only for this clustered case we can 
get maximal MIMO diversity, any other placement of relays 
results in lower diversity gains. Motivated by this fact, we 
will also study the effect of clustering on the relaying systems 
listed above. 

A. Related Work 

Most of the literature on cooperative communications con- 
sider single antenna terminals. The DMT of relay systems 
were first studied in [6] and [39] for half-duplex relays. 
Amplify-and-forward (AF) and decode-and-forward (DF) are 
two of the protocols suggested in [6] for a single relay system 
with single antenna nodes. In both protocols, the relay listens 
to the source during the first half of the frame, and transmits 
during the second half, while the source remains silent. To 
overcome the losses of strict time division between the source 
and the relay, [6] offers incremental relaying, in which there 
is a 1-bit feedback from the destination to both the source and 
the relay, and the relay is used only when needed, i.e. only if 
the destination cannot decode the source during the first half 
of the frame. In [27], the authors do not assume feedback, 
but to improve the AF and DF schemes of [6] they allow the 
source to transmit simultaneously with the relay. This idea is 
also used in [39] to study the non-orthogonal amplify-and- 
forward (NAF) protocol in terms of DMT. Later on, a slotted 
AF scheme is proposed in [40], which outperforms the NAF 
scheme of [39] in terms of DMT. Azarian et al. also propose 
the dynamic decode-and-forward (DDF) protocol in [39]. In 
DDF the relay listens to the source until it is able to decode 
reliably. When this happens, the relay re-encodes the source 
message and sends it in the remaining portion of the frame. 
The authors find that DDF is optimal for low multiplexing 
gains but it is suboptimal when the multiplexing gain is large. 
This is because at high multiplexing gains, the relay needs to 
listen to the source longer and does not have enough time left 
to transmit the high rate source information. This is not an 
issue when the multiplexing gain is small as the relay usually 
understands the source message at an earlier time instant and 
has enough time to transmit. 



MIMO relay channels are studied in terms of ergodic 
capacity in [41] and in terms of DMT in [42]. The latter 
considers the NAF protocol only, presents a lower bound on 
the DMT performance and designs space-time block codes. 
This lower bound is not tight in general and is valid only 
if the number of relay antennas is less than or equal to the 
number of source antennas. 

The multiple-access relay channel (MARC) is introduced 
in [43], [20], [44]. In MARC, the relay helps multiple sources 
simultaneously to reach a common destination. The DMT for 
the half-duplex MARC with single antenna nodes is studied 
in [45], [46], [47]. In [45], the authors find that DDF is DMT 
optimal for low multiplexing gains; however, this protocol re- 
mains to be suboptimal for high multiplexing gains analogous 
to the single-source relay channel. This region, where DDF 
is suboptimal, is achieved by the multiple access amplify and 
forward (MAF) protocol [46], [47]. 

When multiple single antenna relays are present, the pa- 
pers [9], [12], [15], [23], [27], [33], [34] show that di- 
versity gains similar to multi-input single-output (MISO) or 
single-input multi-output (SIMO) systems are achievable for 
Rayleigh fading channels. Similarly, [6], [39], [48], [49] upper 
bound the system behavior by MISO or SIMO DMT if all 
links have Rayleigh fading. In other words, relay systems 
behave similar to either transmit or receive antenna arrays. 
Problem 4 is first analyzed in [50] in terms of achievable rates 
only, where the authors compare a two-source two-destination 
cooperative system with a 2 x 2 MIMO and show that the 
former is multiplexing gain limited by 1, whereas the latter 
has maximum multiplexing gain of 2. 

B. Contributions 

In the light of the related work described in Section II-AI 
we can summarize our contributions as follows: 

> We study Problem 1 with full-duplex relays and compare 
DF and compress-and-forward (CF) [8], [20] strategies 
in terms of DMT for both clustered and non-clustered 
systems. We find that there is a fundamental difference 
between these two schemes. The CF strategy is DMT 
optimal for any number of antennas at the source, the 
destination or the relay, whereas DF is not. 

> We also study Problem 1 with half-duplex relays. This 
study reveals that for half-duplex systems we can find 
tighter upper bounds than the full-duplex DMT upper 
bounds. Moreover, we show that the CF protocol achieves 
this half-duplex DMT bound for any number of antennas 
at the nodes. This is the first known result on DMT 
achieving half-duplex relaying protocols. 

> For Problem 2 we show that the CF protocol achieves a 
significant portion of the half-duplex DMT upper bound 
for high multiplexing gains. Our results for single antenna 
MARC easily extend to multiple antenna terminals. 

• We examine Problem 3 and Problem 4 and develop the 
DMT analysis to understand if the network provides any 
MIMO benefits. Our analysis shows that even for clus- 
tered systems with full-duplex relays, all relay systems 
fall short of MIMO, mainly due to multiplexing gain 
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limitations. The same problem persists in cooperative 
systems with multiple source destination pairs. 

Overall, our work sheds light onto high SNR behavior of 
cooperative networks as described by the DMT, and suggests 
optimal transmission and relaying strategies. 

The paper is organized as follows. Section describes 
the general system model. In Section [nil we give some 
preliminary information that will be used frequently in the 
rest of the paper In Section |IV] we solve the single user, 
single relay problem with multiple antennas for full-duplex 
relays, and in Section |V] we solve the same problem for half- 
duplex relays {Problem 1). Section IVTl introduces MARC, and 
suggests an achievable DMT {Problem 2). In Section IVIII 
we study two problems: the two relay system with a single 
source destination pair {Problem 3), and the two source two 
destination problem {Problem 4). Finally, in Section IVIIII we 
conclude. 

II. General System Model 

For the most general model all the channels in the system 
have independent, slow, frequency non-selective, Rician fad- 
ing. For Rician fading channels, the channel gain matrix is 
written as 



H 



K 



-H,+ 



1 



-H, 



K+1 \l K+l 
where K > 0, Hi and Hg denote the Rician factor, the line 
of sight component and the scattered component respectively. 
The DMT for Rician channels are studied in detail in [51]. In 
[51] the authors find that for finite Rician factor K, the channel 
mean does not affect the DMT behavior, and the system DMT 
will be equal to that of a Rayleigh fading channel with K = 0. 
On the other hand in [52] the authors study the effect of K 
on MISO and SIMO DMT when K approaches infinity. They 
find that for large K, the system diversity increases linearly 
with K. Moreover, when K tends to infinity, the diversity gain 
is infinity for all multiplexing gains up to rank{lli). 

Based on the above observations, without loss of generality, 
in this work we assume a discrete approximation to the Rician 
model: If two nodes are apart more than a threshold distance 
A*, the line of sight component is too weak and the Rician 
factor K can be assumed to be equal to zero. Thus the channel 
gain matrix is distributed as Rayleigh, and we say that the 
nodes are in Rayleigh zones, Fig.fTfa). On the other hand, if the 
inter-node distance is less than A*, the line of sight component 
in the received signal is strong; K can be assumed to be 
infinity and the Rician distribution approximates a Gaussian. In 
this case we say that the nodes are in AWGN zones. Fig. [Tfb). 

For the Rayleigh zone, the channel gain matrix for MIMO 
terminals has independent, identically distributed (i.i.d.) zero 
mean complex Gaussian entries with real and imaginary parts 
each with variance a^. The variance is proportional to 
1/A", where A denotes the internode distance, and a is the 
path loss exponent. If nodes i and j are in the AWGN zone, 
the channel gain matrix from node i to j has deterministic 



entries, all equal to 



and the channel gain matrix has 



. 

\ A 



"\ A 



(M. 



rank 1. There is also a dead zone around the nodes, which 
limits the channel gain. 



(a) (b) 
Fig. 1. Two nodes A and B are in (a) Rayleigli zones, (b) AWGN zones. 



Depending on the locations of the nodes, the Rayleigh or 
AWGN zone assumption results in two important configura- 
tions we will consider: clustered and non-clustered. For the 
clustered system, all the source(s) and some of the relay(s) 
are in the same AWGN zone, and the destination(s) and 
the remaining relay(s) are in another AWGN zone, but the 
source cluster and the destination cluster, which are more than 
the threshold distance A* apart, are in their Rayleigh zones. 
However, for the non-clustered system, every pair of nodes in 
the system are in their Rayleigh zonefl We do not explicitly 
study the systems in which some nodes are clustered and some 
are not in this paper, although our results can easily be applied 
to these cases as well. 

The relay(s) can be full-duplex, that is they can transmit 
and receive at the same time in the same band (Sections IIVI 
and IVIII ). or half-duplex (Sections [V] and [VTll. The transmitters 
(source(s) and relay(s)) in the systems under consideration 
have individual power constraints Pi. All the noise vectors at 
the receivers (relay(s) and destination(s)) have i.i.d. complex 
Gaussian entries with zero mean and variance 1 . Without loss 
of generality we assume the transmit power levels are such 
that the average received signal powers at the destination(s) are 
similar, and we define SNR as the common average received 
signal to noise ratio (except for constant multiplicative factors) 
at the destination. Because of this assumption, for the clustered 
systems we study in Section IVIII the nodes in the source 
cluster hear the transmitters in their cluster much stronger than 
the transmitters in the destination cluster, and for all practical 
purposes we can ignore the links from the destination cluster 
to the nodes in the source cluster. This assumption is the same 
as the level set approach of [11]. For non-clustered systems 
each node can hear all others. 

All the receivers have channel state information (CSI) about 
their incoming fading level^. Furthermore, the relays that 
perform CF have CSI about all the channels in the system. 
This can happen at a negligible cost by proper feedback. We 
will explain why we need this information when we discuss 
the CF protocol in detail in Section |IV] The source(s) does not 
have instantaneous CSI. We also assume the system is delay- 
limited and requires constant-rate transmission. We note that 
under this assumption, information outage probability is still 
well-defined and DMT is a relevant performance metric [53]. 

'Note that all mutual information expressions in the paper will be con- 
sidered as random quantities. However, if two nodes are clustered, the 
channel gains in between these two nodes' antennas take certain values with 
probability one. 

^Because of this assumption, all mutual information expressions in the 
paper should be interpreted as conditioned on the receiver side CSI. We omit 
this conditioning in the expressions for notational simplicity. 



4 



IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. X, NO. X, DECEMBER 2007 



There is also short-term average power constraint that the 
transmitters have to satisfy for each codeword transmitted. For 
more information about the effect of CSI at the transmitter(s) 
and variable rate transmission on DMT we refer the reader 
to [54]. 

III. Preliminaries 

In this section we first introduce the notation, and present 
some results we will use frequently in the paper 

For notational simplicity we write /(SNR)=SNR'^, if 

log/(SNR) 



lim 

SNR^c 



logSNR 



The inequalities > and < are defined similarly. In the rest 
of the paper Ij denotes the identity matrix of size i x i, f 
denotes conjugate transpose, and |.| denotes the determinant 
operation. To clarify the variables, we would like to note that 
Ri denotes the ith relay, whereas i?^ ' denotes transmission 
rates; e.g. i?^^^ will be used for target data rate. 

Let i?^^^ (SNR) denote the transmission rate of the system 
and Pe(SNR) denote the probability of error. Then we define 
multiplexing gain r and corresponding diversity d{r) as 

i?(^)(SNR) 



lim 

SNR^oo log SNR 
log Pe (SNR) 

lim 

SNR^ 



-d(r) 



logSNR 

The DMT of an m x n MIMO is given by d„i„(r), the 
best achievable diversity, which is a piecewise-linear function 
connecting the points (fc,dm„(fc)), where dmn{k) = (m — 
k){n — k), k = 0, 1, min{m, n} [37]. Note that dmn{r) = 

In [37], the authors prove that the probability of error is 
dominated by the probability of outage. Therefore, in the rest 
of the paper we will consider outage probabilities only. 

We know that for any random channel matrix H of size 
n X m and for any input covariance matrix Q of size m x m 
[37], 

sup log |l„ + HQHI'I 

Q>0,traco{Q}<mSNR 

< l0g|l„ +TOSNRHHt| . (1) 

Combined with the fact that a constant scaling in the transmit 
power levels do not change the DMT [37], this bound will be 
useful to establish DMT results. 

In a general multi-terminal network, node k sends informa- 
tion to node I at rate 

V 

where is the number of channel uses, Wki denotes the 
message for node I at node k, and Wki is Wki's estimate 
at node /. Then the maximum rate of information flow from a 
group of sources to a group of sinks is limited by the minimum 
cut [55, Theorem 14.10.1] and we cite this result below. 

Proposition 1: Consider communication among m nodes in 
a network. Let Ci C {1,2, ...,to} and Cf be the complement 
of Cj in the set {l,2,...,m}. Also X*'^*) and X*'^*) denote 



transmitted signals from the sets Ci and Cf respectively. Y'^''* ^ 
denotes the signals received in the set Cf. For information rates 
jliki) fj-om node k to /, there exists some joint probability 
distribution X2, Xm), such that 



E 



for all Ci C {1,2,..., to}. Thus the total rate of flow of 
information across cut-sets is bounded by the conditional 
mutual information across that cut-set. 

We can use the above proposition to find DMT upper bounds. 
Suppose R^'^'''^ ~ r^'^') logSNR denotes the tai-get data rate 
from node k to node I, and r''"'' is its multiplexing gain, 
_R(-^'^i^ — X^fcGC lec ^'"'^'''^ denotes the sum target data rate 



Jki) 



IS Its sum 



across cut-set Ci and r^'-''' = X^fceC i&C 
multiplexing gain. We say the link from fc to Z is in outage if 
the event 

8,1 = {i?(«) < 
occurs. Furthermore, the network outage event is defined as 



N 



U ^ki, 

kJe{l,2,...,m},k^l 



which means the network is in outage if any link is in outage. 
Minimum network outage probability is the minimum value 
of P{£n) over all coding schemes for the network. We name 
the SNR exponent of the minimum network outage probability 
as maximum network diversity, d{f), where f is a vector of 
all r^'^'^'s. Then we have the following lemma, which says 
that the maximum network diversity is upper bounded by the 
minimum diversity over any cut. 

Lemma 1: For each Ci C {1, 2, to}, define the maximum 
diversity order for that cut-set as 



dc,(?-*'^'') = - lim 



logminp(j,^^j 



SNR^oo log SNR 

Then the maximum network diversity d{f) is upper bounded 
as 

d{f) < min{dc.(r('^'))}. 
Proof: We provide the proof in Appendix U ■ 
In addition to Lemma [T] the following two results will also 
be useful for some of the proofs. 

Lemma 2 ([56]): For two n x n positive definite matrices 
A and B, if A — B is positive semi-definite, then \A\ > \B\. 

Lemma 3: For two real numbers x,y > 0, xy/{x + y) < c, 
where c is a non-negative real number, implies x < 2c, or 
y < 2c. Therefore, for two non-negative random variables X 

and Y, P {XY/{X + Y) <c) < P{X < 2c) + P{Y < 2c). 

The proof follows from simple arithmetic operations, which 
we omit here. 

IV. Multiple Antenna Nodes, Single Full-Duplex 
Relay 

The general multiple antenna, multiple source, destination, 
relay network includes the multiple antenna relay channel 
consisting of a single source, destination and relay, as a special 
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Fig. 2. System model for Problem 1, the single relay system. The source, 
the destination and the relay have m, n and k antennas respectively. 



case. Any attempt to understand the most general network 
requires us to investigate the muhiple antenna relay channel 
in more detail. Therefore, in this section we study Problem 
1, in which the source, the destination and the relay has m, 
n and k antennas respectively. This is shown in Fig. |2] As 
clustering has a significant effect on the DMT performance of 
the network, we will look into the non-clustered and clustered 
cases and examine the DF and CF protocols. In this section 
the relay is full-duplex, whereas in Section |V] the relay will 
be half-duplex. 



A. Non-Clustered 

Denoting the source and relay transmitted signals as X5 and 
X/;, when the system is non-clustered, the received signals at 
the relay and at the destination are 



Yr = 
Yd = 



HsdXs 



Zr 

Urd^r 



-'D, 



(2) 
(3) 



where Zr and Zd are the independent complex Gaussian 
noise vectors at the corresponding node. Hsr, Hsd and Urd 
are the kxm, nxm and nx k channel gain matrices between 
the source and the relay, the source and the destination, and 
the relay and the destination respectively. 

Theorem 1: The optimal DMT for the non-clustered system 
of Fig. 121 d{r), is equal to 

d{r) = min{d„j(„+fe)(r),d 

and the CF protocol achieves this optimal DMT for any m, n 
and k. 

Proof: 

1 ) Upper Bound: The instantaneous cut-set mutual information 
expressions for cut-sets Cs and Cd are 



Ics = I{Xs;YrYd\Xr) 



(4) 
(5) 



To maximize these mutual information expressions we need 
to choose Xs and X/f complex Gaussian with zero mean 
and covariance matrices having trace constraints P5 and Pr 
respectively, where P5 and Pr denote the average power 
constraints each node has [57]. Moreover, the covariance 
matrix of X5 and X^; should be chosen appropriately to 



maximize Iq^ . Then using ([TJ to upper bound Iq^ with 
and Iq^ with /^^ we can write 



where 



KsR,D 



Ics < lis ^^og K's^j^a 

ICo < I'cn '=^'^°SKsR.D, 



Ifc+n + tls,RDtis,RDPs 



In + JlsR,DnlRDiPs + PR) 



with 

Hs.RD 



H-SD 



H 



SR,D 



SD 



H 



RD 



Cs 

(6) 
(7) 

(8) 
(9) 



(10) 



The above bounds suggest that the CSI at the relay does 
not improve the DMT performance under short term power 
constraint and constant rate operation. The best strategy for the 
relay is to employ beamforming among its antennas. For an m 
antenna MIMO, with total transmit power P, the beamforming 
gain can at most be mP [37], which results in the same DMT 
as using power P. Therefore, CSI at the relay with no power 
allocation over time does not improve the DMT, it has the 
same the DMT when only receiver CSI is present. 

Note that P(/^J=SNR"'*^. ^''^ i = S,D, with 4^.(r) = 
d,n(n+k}i'r) and d'^^ir) = d(„+fe)„(r). Then using Lemma [T] 
one can easily upper bound the system DMT by 

d{r) < miii{d„(„+fe)(r),d(,„+fe)„(r)}, 

for a target data rate R^"^^ — r log SNR. 
2 ) Achievability : To prove the DMT upper bound of Theorem[T] 
is achievable, we assume the relay does full-duplex CF as we 
explain below. We assume the source, and the relay perform 
block Markov superposition coding, and the destination does 
backward decoding [20], [58], [59]. The encoding is carried 
over B blocks, over which the fading remains fixed. In the CF 
protocol the relay performs Wyner-Ziv type compression with 
side information taken as the destination's received signal. For 
this operation the relay needs to know all the channel gains 
in the system. 

For the CF protocol, as suggested in [8] and [20], the relay's 
compression rate has to satisfy 



I{Yr-Yr\^rYd)<IO^r.Yd), 



(11) 



in order to forward reliably to the destination. Here 
denotes the compressed signal at the relay. The destination can 
recover the source message reliably if the transmission rate 
— ^ \Qg SNR of the source is less than the instantaneous 
mutual information R^^-^^ 



=/(X5;Y^Yz,|Xfl). 



(12) 



We assume X5 and 'Kr are chosen independently, and have 
covariance matrices l„iPs/ra and IfeP^/fc respectively. Also 
Yfl = Y r + Zr, where Zr is a length k vector with complex 
Gaussian random entries with zero mean. Zr has covariance 
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matrix Njflk, and its entries are independent from all other 
random variables. We define 

m 



Ls,D = 



HrdH 



m 



RD"-R.D- 



R 



k 



RD — 



■ TO 
' TO 



(iVij + l)Ife 






In 



^SMD — 

Then we have 

I{YR;Yn\Xj,YD) = 

I{Xn;YD) - 



I 



/c+n 

log 
log 



(13) 
(14) 
:-5) 
(16) 



L 



S.RD 



LsB.,D 



L 



S.D 



To satisfy the compression rate constraint in (fTTT i. using the 
CSI available to it, the relay ensures that the compression noise 
variance Nji satisfies Nji = \/ Ls.rd/ Lsr.d- Note that both 
sides of this equation are functions of Nb.- Then 

- /(Xs;Y«Yc|X«) 

Ls,RD 



= log- 



= log 



1 



\/Ls,RD {/LsR,D 



(17) 



\ {/Ls.RD + yLsRJD ) 

To prove the DMT of ( [TtI i we need to find how probability 
of error decays with increasing SNR when the target rate 
increases as i?^-'"^ = r log SNR. As the error events are 
dominated by outage events, we use the following bound on 
the probability of outage 



P (outage at D) 

= P f/(Xs; YfiYz,|X;^) < rlogSNR) 



P 



\/Ls,RD {/LSR,D 

{/Ls.RD + {/L 



< SNR- 



SR.D 



(a) 
< 



(b) 
< 



P 



^'s.RD 



SR,D 



< SNR^^ 



P 



SMD + yLsRJ) 

s,RD < 2SNR^ 

+ p(Vp^<2SNR* 

P {L's^nn < 2'^SNR'-) 
+ P (LsR.D < 2*'-SNR'-) 



(18) 
(19) 

(20) 

(21) 

(22) 
(23) 



= SNR- 



SNR- 



where for (a) we first used Lemma |2] to show Ls.rd > 
^'s RD ^^'^ '■1'^ ^^^^ '■l'^'- '■1'^ islio xy/{x + y) is monotone 
decreasing with decreasing x for x, j/ > 0, (&) follows from 
Lemma |3] and (c) follows because logL'g and \ogLsR^D 



are same as the cut-set mutual information expressions /cs 
and except a constant scaling factor of SNR, and a 
constant scaling in SNR in the probability expression does 
not change the diversity gain. We conclude that the system 
DMT dcF{r) > min{d„(„+fe) (r), d(„+s.)„(r)}. This resuh 
when combined with the upper bound results in 



dcF{r) 



T^{dm{n+k){r),d(^m+k)n{r)-} 



d-Dpir) = 



(24) 



As an alternative to the CF protocol, the relay can use the 
DF protocol. When the source, the destination and the relay 
all have a single antenna each, it is easy to show that the DF 
protocol also achieves the DMT upper bound, which is equal 
to di2{r). The following theorem derives the DMT of the DF 
protocol for arbitrary to, n and k and shows that the optimality 
of DF does not necessarily hold for all to, n and k. 

Theorem 2: For the system in Fig.|2] DF achieves the DMT 

dmn (^) ~t~ dmk 

if < r < min{m, n, k} 

if min{TO, n, fc} < r < min{m, n\ 

Proof: We provide the proof in Appendix ■ 

We next consider examples for the DF DMT performance 
and compare with Theorem [1] If to or n (or both) is equal 
to 1, we find that DF meets the bound in Theorem [T] and is 
optimal irrespective of the value of fc. Similarly we can show 
that for cases such as (to, n, fc) = (3, 2, 2) or (m, n, fc) — 
(4,2,3), as d(,n+k)n{r) < d,nn{r) + d„rk{r) for all r, DF is 
optimal. A general necessary condition for DF to be optimal 
for all multiplexing gains is to > n. If to < n, then dmn (r) + 
d,nk{r) < dm(n+k){r) < d(^m,+k)n, and DF will be suboptimal. 

Whenever minjm, n, fc} = fc, the degrees of freedom in the 
direct link is larger than the degrees of freedom in the source to 
relay link, that is min{TO, n} > min{TO, fc}. For multiplexing 
gains in the range min{TO, n, k} < r < min{TO, n}, the relay 
can never help and the system has the direct link DMT dmn(r)- 
Therefore, DF loses its optimality. For example, if (m, n, fc) ~ 
(3, 2, 1), then DF is optimal only for multiplexing gains up to 
1/2, but for 1/2 < r < 2, DF is suboptimal. In particular, DF 
does not improve upon ^32(7') in the range 1 < r < 2. 

Fig. [3] shows the CF and DF DMT for (to, n, fc) = (2, 2, 1), 
and Fig. |4] shows the CF and DF DMT for (to, n, fc) = 
(2,2,2). When we compare the figures, we see that the CF 
protocol is always DMT optimal, but the DF protocol can 
still be suboptimal even when the source to relay link has the 
same degrees of freedom as the link from the source to the 
destination. The suboptimal behavior of DF arises because the 
outage event when the relay cannot decode can dominate for 
general to, n and fc. In addition to this, for multiplexing gains 
larger than min{TO, n, fc}, the relay never participates in the 
communication because it is degrees of freedom limited and 
cannot decode large multiplexing gain signals. For this region, 
we observe the direct link behavior. We conclude that soft 
information transmission, as in the CF protocol, is necessary 
at the relay not to lose diversity or multiplexing gains. 

Fig. |5] shows the outage probability versus total SNR for 
DF and CF protocols for (m, n, fc) = (2,2,1), R^'^'> = 
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-B- Tradeoff for CF (Optimal) 
-e- Tradeoff for DF 







0' • ^ 

1 2 



r 

Fig. 4. The source has 2, the destination has 2, and the relay has 2 antennas, 
(m, n, k) = (2, 2, 2). The network is non-clustered. 

rlog{Ps + Pr), r — 1.5. The channel gain matrices Hsr^, 
HsD and Urd have i.i.d. complex Gaussian entries, with 
real and imaginary parts zero mean and variance 1/2 each. 
We have P5 = 2Pr. The figure also includes the 3x2 
MIMO for comparison. We assume the total power constraint 
is the same for both the MIMO and relay systems. In the 
MIMO system the antennas share the total power equally and 
send uncorrected signals. We observe that while DF achieves 
d = 0.5, CF achieves d = 1 and performs similar to 3 x 2 
MIMO as predicted by Theorem [T] 

The above analysis also reveals that CF and DF protocols 
do not always behave similar, unlike the single antenna relay 
system. The degrees of freedom available also has an effect 
on relaying strategies. 

B. Clustered 

Clustering can sometimes improve the system performance, 
since it eliminates fading between some of the users. We 
will observe an example of this in Section IVIII when there 



are multiple relays. Therefore, in this subsection we study 
the DMT behavior of a single relay system, when the relay 
is clustered with the source. The analysis presented in this 
subsection can easily be modified if relay is clustered with 
the destination. The system input and output signals are same 
as ^ and (O but for the clustered case all the entries of Hsr 
are equal to \/G. 

Theorem 3: For the system in Fig. |2] when the relay is 
clustered with the source, the CF protocol is optimal from 
the DMT perspective for all (to, n, k). 

We omit the proof, as the achievability follows the same 
lines as in Theorem [T] and results in the same outage proba- 
bility expression in ( |23] |, which is equal to the upper bound. 

We next compute the DMT of the clustered system explicitly 
for m = 1, 2, arbitrary n and k. We conjecture the same form 
holds for arbitrary to as well. 

Theorem 4: For the clustered system of Fig. |2] for 
(to, n, k), 171—1,2, the DMT is given by 

{d{m+k)n{r) ifO<r<l 
min{dm(„+i)(r), d(™+fc)„(r)} 
if 1 < r < min{TO, n} 

Proof: We provide the proof in Appendix |lll] ■ 
For the clustered case, for any to or k, Hsr has rank 1, hence 
we have the following conjecture. 

Conjecture 1: Theorem |4] is true for arbitrary {m,n,k). 

We observe that if m = 2, although the source and 
the relay both have multiple antennas, as the channel gain 
matrix in between is AWGN and has rank 1, it can only 
support multiplexing gains up to 1. This is because having 
multiple antennas at the transmitter and/or the receiver in 
an AWGN channel only introduces power gain. Therefore, 
Ics, the mutual information across cut-set Cs, never results 
in outage for multiplexing gains up to 1. For multiplexing 
gains r > 1, this cut-set results in a DMT of (i2(n+i)(^)' 
even though the relay has k antennas. The next theorem is a 
counterpart of Theorem |2] for the clustered case. 
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-Tradeoff for CF (Optimal) 
■Tradeoff for DF 



Comparing with the upper bound, we can see that clustering 
improves the DF performance in the range < r < 1, where 
DF achieves the upper bound. However, for muhiplexing gains 
larger than 1, DF is still suboptimal. In fact, in this range 
the relay can never decode the source even though they are 
clustered and hence cannot improve the direct link perfor- 
mance. Although clustering improves the DF performance 
for low multiplexing gains, it is not beneficial for multiple 
antenna scenarios in terms of DMT, it can in fact decrease the 
optimal diversity gain. This is because when two nodes have 
multiple antennas, clustering decreases the degrees of freedom 
in between. This can also be observed comparing Theorem H] 
and Theorem[Tl as well as the optimal strategies in Fig.|2]with 
Fig. 13] We will also study the effects of clustering in single 
antenna multiple relay scenarios in Section IVIII 



Fig. 6. The source has 2, the destination has 2, and the relay has 1 antenna, 
(m, 71, k) = (2, 2, 1). The relay is clustered with the source. 



-B- Tradeoff for CF (Optimal) 
X. -e- Tradeoff for DF 




^^^^^^ 



(r) 



Fig. 7. The source has 2, the destination has 2, and the relay has 2 antennas, 
(m, n, k) = (2, 2, 2). The relay is clustered with the source. 



Theorem 5: For the system in Fig. |2l when the relay is 
clustered with the source, the DF protocol achieves the DMT 

d{,n+k)n{r) if < r < 1 
dmnir) 1 < r < min{m,n} 

Proof: The outage probability for DF is the same as ( |4TI ). 
in the non-clustered case of Appendix Ull If < r < 1, then 
the probability that the relay is in outage is 0. On the other 
hand, if 1 < r < min{m, n}, the probability that the relay can 
decode is 0, since the source-relay channel can only support 
multiplexing gains up to 1. ■ 

We have seen in Section IIV-AI that DF is in general subop- 
timal for non-clustered multi-antenna relay channel. However, 
once we cluster the relay with the source, there are no more 
outages in the source-relay channel for multiplexing gains up 
to 1, and the DF performance improves in this range. However, 
even with clustering DF does not necessarily meet the DMT 
upper bound for arbitrary (m, n, k). 

Fig. |6] compares the clustered CF and DF DMT for 
{m,n,k) = (2,2,1), and Fig. [T] for {m,n,k) = (2,2,2). 



V. Multiple Antenna Nodes, Single Half-Duplex 
Relay 

In the previous section, we studied the relay channel when 
the relay is full-duplex. Although this is an ideal assumption 
about the relay's physical capabilities, it helps us understand 
the fundamental differences between the DF and CF protocols. 
In this section we assume a half-duplex, non-clustered relay to 
study how this affects the DMT behavior of the relay channel. 

In half-duplex operation a state variable Q, which takes 
the value qi if the relay is listening, or q2 if the relay is 
transmitting, controls the relay operation. For a more general 
treatment that considers three different states depending on 
whether the relay is in sleep, listen or talk states see [60]. Our 
results in this section would also be applicable for this case 
as well. 

Depending on how the state Q is designed, half-duplex 
protocols can be random or fixed. In fixed protocols, the state 
does not convey additional information to the destination via 
the state random variable Q, whereas in random protocols the 
relay breaks its transmission and reception intervals into small 
blocks to send extra information through the state. This is 
equivalent to considering the random binary state as a channel 
input and designing code books to convey information through 
Q. 

Another categorization based on the state variable Q is 
dynamic versus static. If the state is controlled based on 
channel realizations, we have a dynamic protocol. On the other 
hand, if Q does not depend on CSI, the protocol is called 
static. Note that fixed protocols are included in random ones, 
and static protocols in dynamic ones. The most commonly 
used relaying protocols are fixed and static, and of the form 
shown in Fig. [8] The DDF protocol of [39] is an example to 
a fixed, dynamic protocol. 

For the multiple antenna half-duplex relay channel, using 
Lemma [T] directly, results in the full-duplex bound, which is 
not tight for half-duplex operation. Therefore, we first state the 
following lemma to provide a half-duplex DMT upper bound 
for random, static protocols. The lemma also suggests that 
sending information through the state does not improve DMT. 
LemmalUcan be modified for random, dynamic state protocols 
as well. 
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Fig. 8. A fixed and static lialf-duplex relaying protocol, where the relay 
listens for t fraction of time, transmits for (1 — t) fraction, where t is a 
constant. 



Lemma 4: For the multiple antenna half-duplex relay chan- 
nel, the half-duplex DMT upper bound for random, static state 
protocols is equal to 



d{r) < mayi min{dcs{r,p{q)),dco{r,piq))}, 
pyi) 



(25) 



where 



dCi{r,p{q)) 



= — lim 

SNR^oo 



logSNR 



(26) 



i^S,D. 



Proof: We provide the proof in Appendix |IVl ■ 
Our next theorem and corollary provide the first half-duplex 
DMT achieving relaying protocol in the literature. 

Theorem 6: For the random, dynamic state, half-duplex 
relay channel with m antenna source, k antenna relay and 
n antenna destination, the CF protocol is DMT optimal. 

Corollary 1: For {m,n,k) — (1,1,1), the half-duplex 
DMT upper bound is equal to the full-duplex DMT, 2(1 — r). 
Therefore, CF is a DMT optimal half-duplex protocol for the 
single antenna relay channel. 

Proof: [Theorem First, we prove that CF is optimal 
among static protocols and then show that the same proof 
follows for dynamic protocols as well. 

At state qi, the received signals at the relay and the 
destination are 

Y_R,i — lIsB.^s,i + 
Yd,i = Hs_dXs,i + Zd,i 

and at state q2, the received signal at the destination is given 
as 



■D,2 



Here X5 ;, and X^^ ; are of size m, and k column vectors 
respectively and denote transmitted signal vectors at node S, 
and R at state qi, 1 — 1,2. Similarly Y/j ; and ; are the 
received signal vectors of size k and n. 

We first find an upper bound to the DMT using Lemma |4] 
Without loss of generality we use a fixed state static protocol 
as shown in Fig. [8] This is justified by the proof of Lemma H) 
which states that fixed and random protocols have the same 
DMT upper bound. For the half-duplex relay channel using the 
cut-set Cs around the source and Cd around the destination 



as shown in Fig. |2] we have [19] 

Icsit) - tI{Xs;YnYD\qi) 

+ (l-t)/(X5;Yz, 1X^^,92) 

Ic^it) = tI{Xs;YD\qi) 

-f (l-t)/(X5Xfl;Yc|g2). 

We define 



K 



S,D 



(27) 
(28) 

(29) 



Then we can upper bound Ics (t) and /c^ {t) with Iq^ (t) and 

Ics (t) < I'cs (*) = t log K's,RD + (1 - t) log Ks.D (30) 

ICo (t) < I'cn (*) = t log Ks,D + (1 - t) log KsR,D (31) 

where Kg and Ksr.d are defined in ^ and Q. 

For a target data rate i?^-'") = r log SNR, and for a 
fixed t, if P(/c^(<) < i?(^))=SNR"'^^:.(''^*\ i ^ S,D, then 
dci{r,t) of Lemma H] satisfies dci{r,t) < d'^Xr^t), where we 
denoted dci{r,p{q)) with dci{r,t) with an abuse of notation. 
Therefore, the best achievable diversity for the half-duplex 
relay channel for fixed t satisfies 



d{r,t) < mm{d'c^{r,t),d'cjr,t)}. 



(32) 



Optimizing over t we find an upper bound on the static 
multiple antenna half-duplex relay channel DMT as 



d{r) < maxmin{4s {r, t),d'c^ {r, t)}. 



(33) 



Appendix [V] shows that half-duplex CF achieves the upper 
bound in (|33] ). For dynamic protocols, the DMT upper bound 
will change because of the CSI available at the relay. Ap- 
pendix |V] also shows that if CF is allowed dynamic operation, 
it achieves the dynamic DMT upper bound as well. ■ 

A. Static Half-Duplex DMT Computation 

In general it is hard to compute the exact DMT of The- 
orem |6] In particular for static protocols, to find (r, t) 
and (r, t) for general to, n and k we need to calculate 
the joint eigenvalue distribution of two correlated Hermitian 
matrices, H^dH^^ and Hs.Roiig jijj or HgcH^^ and 
Hsjijjiigj^ jj. However, when to = 1, both Hg^iH^^ 



t 

S,RD 



reduce to vectors and it becomes easier 



and Hs RD^ 

to find d'(,^{r,t). Similarly, when n — 1, HsdH^^ and 
Hsfl dH^^ are vectors, and d'^^ (r, t) can be found. 

An explicit form for (r, t) is given in the following 
theorem. 

Theorem 7: For m — 1, d[,^ (r, t) is given as 



n + k — kj 
{n + fc)(l - r) 



if r < t, md t < ^ 

— ' — n-\-k 

if r > t, and i < 



if t > 



k 

n+k 



For n = 1 and for arbitrary to and k, d'^^ (r, t) has the same 
expression as d'^^ (r, t) if n and t are replaced with to and 
(1 — i) in the above expressions. 
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Proof: The proof follows the ideas presented in [37], 
[39], and is provided in Appendix I VII ■ 

Although we do not have an explicit expression for (r, t) 
or d'^^ {r, t) for general (m, n, k), we can comment on some 
special cases and get insights about multiple antenna, half- 
duplex behavior First we observe that d'^^ {r, t) and (r, t) 
depend on the choice of i, and the upper bound of ( [33] l is not 
always equal to the full-duplex bound. As an example consider 
{■m,n,k) = (1,1,2), for which d'i^^{r,t) is shown in Fig. |9] 
To achieve the full-duplex bound for all r, d'^^ {r, t) needs 
to have t > 2/3, whereas d'Q^{r,t) needs t < 1/3. As both 
cannot be satisfied simultaneously, d{r, t) will be less than the 
full-duplex bound for all t. 

On the other hand, to maximize the half-duplex DMT it is 
optimal to choose t — \/2 whenever m — n.To see this, we 
compare ( [30l l with (|3TT i, and note that both K'g j^jj > Ks,d 
and KsR.D > Ks,d for m = n. Furthermore, for m — n 
d'i^^{r,t) = dQ^{r,l — t), and d'('^{r,t) is a non-decreasing 
function in t. Therefore mm{d'i^^{r,t),dQ^{r,t)} must reach 
its maximum at t = 1/2. 







\ >. 
I >^ 


1>= 2/3 

1< = 2/3 


\ N. 
\ 

\ >v 

\ N. 

\ ^ 
~ ~ ^ 





t r 1 



Fig. 9. DMT upper bound for the cut-set around the source, Cg . The source 
has 1, the destination has 1, and the relay has 2 antennas, (m, n, fc) = 
(1,1,2). The network is non-clustered. Note that as m = n, d'^^(r,t) = 
d'^^ (r, 1 — t). The upper bound in 133) reaches its maximum for t = 1/2. 
The sohd line in the figure is also equal to the full-duplex bound. 



B. Discussion 

When {m,n,k) = (1,1,1), the best known half-duplex 
DMT in the literature is provided by the dynamic decode- 
and-forward (DDF) protocol [39]. The DDF protocol achieves 

/ 2(1 - r) if < r < i 
dDDF[r) ==<^l-r ifi<r<l' 

which does not meet the upper bound for \ < r < \, 
as in this range, the relay does not have enough time to 
transmit the high rate information it received. We would like 
to note that, if the relay had all CSI, the DMT of the DDF 
protocol would not improve. With this CSI the relay could 
at best perform beamforming with the source; however, this 
only brings power gain, which does not improve DMT. It is 
also worth mentioning that when only relay CSI is present, 
incremental DF [6] would not improve the DMT performance 
of DF. Unless the source knows whether the destination has 
received its message or not, it will never be able to transmit 
new information to increase multiplexing gains in incremental 
relaying. 

In general it is hard to compute the DMT of multiple 
antenna DDF. This is because the instantaneous mutual in- 
formation DDF achieves in a multiple antenna relay channel 
is equal to Ico (t) of (|28] | where t is the random time instant 
at which the relay does successful decoding. Thus it is even 
harder to compute the DMT for this case than for fixed t. 
Moreover, we think that the multiple antenna DDF perfor- 
mance will still be suboptimal. In Section |IV] we showed that 
for a multiple antenna full-duplex relay system, the probability 
that the relay cannot decode is dominant and the DF protocol 
becomes suboptimal. Therefore, we do not expect any relay 
decoding based protocol to achieve the DMT upper bound in 
the multiple antenna half-duplex system either. This conjecture 
is also demonstrated in Fig. [TO] which shows the outage 
probability versus total SNR for DDF and CF protocols for 
(m,n,fc) = (2,2,1), = rlog{Ps + Pr), r = 1.5, 




10 20 30 40 50 

total SNR (dB) 



Fig. 10. Outage probability versus total SNR for the non-clustered multiple- 
antenna, single half-duplex relay system, {m,n,k) = (2,2,1), r = 1.5, 
t = 0.5. 

t = 0.5. Source has twice the power relay has. The matrices 
VLsR, tlsD and H/^x) have i.i.d. complex Gaussian entries 
with real and imaginary parts zero mean and variance 1/2. 
We observe that the diversity gain the CF protocol achieves is 
approximately 0.90, whereas the DDF protocol approximately 
achieves 0.47. 

VI. The Multiple-Access Relay Channel 

The most general network we introduced in Section I] 
includes the multiple access relay channel (MARC) as a 
subproblem, {Problem 2). The model for MARC is shown 
in Fig. [TT] Our emphasis is on half-duplex MARC. As in 
Section [V] without loss of generality we consider a static, fixed 
state protocol, where the relay listens for t fraction of time, 
transmits for (1 — t) fraction and sources transmit all the time. 
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at state qi (when the relay listens) and at state q2 (when the 
relay transmits), the received signal at the destination is given 
as 



Y 



D.2 — 



H 



SiD^Sia 



H 



H 



■'D.2- 



Here X^j j, Xs^ j and Xj^ ^ are of size mi, m2 and k column 
vectors respectively and denote transmitted signal vectors at 
node Si, S2 and R at state qi, for I = 1,2. Similarly Y^^i 
and Yd, I are the received signal vectors of size k and n. 
UsD, iisiD, iis2D, HsiR, Hg^fl, Urd are the channel gain 
matrices of size n x nii, n x m2, k x nii, k x 7712, n x k 
respectively. The system is non-clustered. 

In this section we examine the DMT for the MARC. We 
present our results for the MARC with single antenna nodes 
to demonstrate the basic idea. 

The DMT upper bound for the symmetric MARC occurs 
when both users operate at the same multiplexing gain r/2, 
f = (r/2, r/2), and is given in [45] as 



dMARcir) < 



3(1 -r) 



if < r < i 

2<r<^ 



if 1 



(34) 



which follows from cut-set upper bounds on the information 
rate. Although this upper bound is a full-duplex DMT bound, 
it is tight enough for the half-duplex case when each node 
has a single antenna. We see that this upper bound has the 
single user DMT for r < ^, and has the relay channel DMT 
with a two-antenna source for high multiplexing gains. This is 
because for low multiplexing gains, the typical outage event 
occurs when only one of the users is in outage, and at high 
multiplexing gains, the typical outage event occurs when both 
users are in outage, similar to multiple antenna multiple-access 
channels [61]. 

In Sections |IV] and |V] we have observed that CF is DMT 
optimal for full-duplex and half-duplex multi-antenna relay 
channels. This motivates us to study the performance of CF 
in MARC. 

^This should not be confused with the notation of [61], in which r denotes 
the per user multiplexing gain in case of symmetric users. 



Theorem 8: For the single antenna, half-duplex MARC the 
CF strategy achieves the DMT 

2(1- 



if < r < I 

dMARC.CF{r) ^ { l-§ if |<r<| 



if 



< r < 1 



3(1 -r) . , 

This DMT dMARC,CF{r) becomes equal to the upper bound 
for r > 4/5. 



Proof: The proof is provided in Appendix IVIII ■ 
To achieve the above DMT performance, two types of 
operation are necessary. For low multiplexing gains, < r < 
2/3, and ^2 utilize time sharing, and equally share the 
relay. Here both Si and ^2 transmit for the half of the total 
time, for 1 /4 of the whole time slot R helps 5*1 only, and in the 
last quarter, R helps 5*2. Then we can directly apply the results 
obtained in Section IVl which results in the DMT 2(1 — r) in 
terms of the sum multiplexing gain. For high multiplexing 
gains, 2/3 < r < 1, both sources transmit simultaneously. In 
this multiple access mode, for 4/5 < r < 1, both users being 
in outage is the dominant outage event, the system becomes 
equivalent to the multiple antenna half-duplex relay channel, 
and CF achieves the DMT upper bound. 

For comparison the achievable DMT with DDF for MARC 
satisfies [45]: 

f2-r if 0<r<i 
dMARC.DDF{r/2,r/2) >\ 3(1 - r) if i < r < § . 

21^ if ^ <r <1 

\ r 3 — — 

We also compare our results with the MAF protocol for the 
MARC channel [46], [47] in Fig. [12] The MAF performance 
is given as 



dMARC,MAFir/2,r/2) = 



2-3 
3(1- 



if 
if 



0<r<i 
|<r<l 



We observe that for low multiplexing gains, when single 
user outage is dominant, it is optimal to decode the sources; 
however for high multiplexing gains, compression works bet- 
ter The MAF protocol is also optimal for high multiplexing 
gains. 
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Fig. 13. System model for Problem 3, the single source-destination two relay 
system, each node has a single antenna. 

In Section |IV] we observed that for a full-duplex relay 
channel, when terminals have multiple antennas, DF becomes 
suboptimal, whereas CF is not. Hence, we conjecture that 
DDF will not be able to sustain is optimality even in the low 
multiplexing gain regime when the terminals have multiple 
antennas. Moreover, it is not easy to extend the MAF protocol 
for multiple antenna MARC. Even when we have one source, 
the DMT for the multiple antenna NAF protocol for the relay 
channel is not known, only a lower bound exists [42]. On 
the other hand, for the multiple antenna case CF will still 
be optimal whenever decoding all sources together is the 
dominant error event. However, for some antenna numbers mi, 
1712, and 71, single-user behavior will always dominate [61]. 

VII. Single Antenna Nodes, Multiple Relays 

In this section we examine Problem 3 and Problem 4 to see 
how closely a cooperative system can mimic MIMO in terms 
of DMT. We first study a single source destination pair with 2 
relays (Problem 3) in Section lVII-AI then consider two sources 
and destinations (Problem 4) in Section [VII-BI In both cases, 
each node has a single antenna. We also assume the nodes are 
full-duplex so that we can observe the fundamental limitations 
a relaying system introduces. 

A. Single Source-Destination, Two Relays 

In this system there is a single source-destination pair and 
two relays as shown in Fig. [T3] The channel is characterized 
by 

YR-^ = asR^Xs + fiR^ji^Xi^,^ + (35) 
= hsR^Xs + Iir-^r^Xr^ + Zr^ (36) 
Yd ~ hsD^s + It-r^dXr^ + ur^dXr., + Zd (37) 

where Xi and Yi, i — S,Ri,R2,D, are transmitted and 
received signals at node i respectively. The channel gains 
hij, i,j = S, Ri, R2, D, are independent, zero mean complex 
Gaussian with variance 2(T^, where cr^ is defined in Section HIl 
As discussed in Section |ll] we assume the R2 to i?i link 
hR2Ri, which is the dashed line in Fig. [13] is present only 
if the system is not clustered. If the system is not clustered, 
then the channel gains are also Rayleigh. On the other 




Fig. 14. DMT for the single source-destination, two relay system, each node 
has a single antenna. 

hand, if the system is clustered, then asR^ and ur^d are 
equal to ^JGsr^ and Gr^d respectively, which are the 
Gaussian channel gains. Z, denotes the AWGN noise, which 
is independent at each receiver The source, the first relay, 
and the second relay, R2, have power constraints P5, 
Pfjj and Pr^ respectively. We assume the target data rate 
jl(T) _ rlogSNR. The following theorems summarize the 
main results of this section. 

Theorem 9: The optimal DMT for the non-clustered system 
of Fig. [13] d{r), is equal to 

d{r) = di3(r). 

This optimal DMT is achieved when both relays employ DF 
strategy. 

Proof: Please refer to Appendices IVIIII and |IX] for the 
DMT upper bound and achievability results respectively. ■ 
Theorem 10: The optimal DMT for the clustered system 
of Fig. [T3I where Ri is clustered with the source and R2 is 
clustered with the destination, d{r) is equal to 

. . U22(r) if r<l 

\ if r > 1 • 

The mixed strategy, where i?i does DF and R2 does CF 
achieves the optimal DMT. 

Proof: Please refer to Appendices IVIIII and |X] for the 
DMT upper bound and achievability results respectively. ■ 

Theorem |9] says that if the system is non-clustered it can at 
most have a transmit or a receive antenna array DMT behavior, 
but cannot act as a MIMO in terms of DMT. On the other hand. 
Theorem [TO] confirms the fact that the multiplexing gain for 
the clustered system is limited by 1. However, for all r < 1, 
the clustered system can mimic a 2 x 2 MIMO, which means 
d = 1 is achievable at r = 1. 

The DMT performances for non-clustered and clustered 
systems as well as 2 x 2 MIMO are illustrated in Fig. [14] 
We also display the outage probability versus SNR for this 
clustered case in Fig. [15] for i?'^^ — r \og{Ps + Pr^ + ^'^2)' 
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Fig. 15. Outage probability for the clustered single source-destination, two 
relay system, each node has a single antenna. Gsr^ = Gji^u = 10, r = 1. 



where r = 1. We assumed GsRi — Gr^d — 10 and 
Ps = = IOPr^. Also, hij, i,j = S,Ri,R2,D, are i.i.d. 
with cr^ ~ 1/2. For comparison, we also show the outage 
probability of a 2 x 2 MIMO channel, where the 2 transmit 
antennas share the total power equally and send uncorrected 
signals. The MIMO channel and the relay system have the 
same total power constraint. We observe that as predicted, the 
clustered relay network has the same diversity as the 2x2 
MIMO and at r = 1, d = 1 is achievable. 

We would like to note that for the clustered case CF is 
essential at i?2, and a strict decoding constraint at i?2 would 
limit the system performance. If both relays do DF, Ri will 
always be able to decode for all multiplexing gains < r < 1, 
as the S — Ri channel can support rates up to \og{l + GsRiP)- 
Thus, it is as if there is a two-antenna transmitter However, 
i?2 may or may not decode. Adapting Appendix |IX] to the 
clustered case we can easily find the probability of outage 
at the destination from (|54]i as P(outage at D)=SNR"'^'^i'^''\ 
which shows that the decoding constraint at R2 limits the 
system performance, the system still operates as a transmit 
antenna array. Even though one could improve upon this 
strategy by using the DF protocol of [20], which allows 
the relays to process the signals they hear from the source 
and the other relay jointly, this still does not provide 2x2 
MIMO behavior In this case both the destination and R2 
observe 2 x 1 DMT, and P(outage at D)=^SNR"2'^=i('") = 
SNR"^^^^*^', which is still suboptimal as it cannot achieve 
the upper bound of Theorem [TO] Although the destination 
can always understand R2 reliably (because of the clustering 
assumption), whenever both of them fail, the system is in 
outage. However, for the receive cluster, CF fits very well. 
If the received signal at the destination has high power due 
to large hsD and hn^n, then R2 to destination channel has 
lower capacity because in the decoding process Yd is treated 
as interference. On the other hand, the correlation between the 
relay and destination signals is higher and a coarse description 
Yfl2 is enough to help the destination. However, if the side 
information has low received power, the R2 to destination 




Fig. 16. System model for Problem 4, the two-source two-destination system, 
each node has a single antenna. 



channel has higher capacity and R2 can send the necessary 
finer information as the correlation is less. 

In [20, Theorem 4], the authors prove an achievable rate for 
a multiple relay system, in which some of the relays DF and 
the rest CF. Furthermore, the relays that perform CF partially 
decode the signals from the relays that perform DF. Performing 
this partial decoding leads to higher achievable rates. However, 
to achieve the DMT upper bound, for both the non-clustered 
and clustered cases, there is no need for partial decoding and 
a simpler strategy is enough. 

Note that same multiplexing limitations in Theorem [10] 
would occur when the source has two antennas and a single 
antenna relay is clustered with a single antenna destination or 
the symmetric case when the destination has two antennas and 
a single antenna relay is clustered with a single antenna source. 
These multiple antenna, single source-destination, single relay 
cases were discussed in detail in Section |IV] In addition to 
these, we investigate whether the multiplexing gain limitation 
is due to the fact that there is only a single source-destination 
pair in the next subsection. 

B. Two-Source Two-Destination Cooperative System 

We consider two sources and two destinations, where 
sources cooperate in transmission and destinations cooperate 
in reception {Problem 4). The system model is shown in 
Fig. [16] Problem 3 studied in Section lVII-Al would be a special 
case of this, if one source has no information to send. 

First we examine the multi-cast scenario, when both desti- 
nations are required to decode both sources. This is analogous 
to MIMO systems and represents the information transfer from 
a group of antennas to another group of antennas. We define 
individual target data rates ' — ri log SNR and R^'^^^ ) = 
r2 log SNR, with a sum target data rate of i?'^' — r log SNR, 
r — ri + r2, r = (ri, r2). Using the cut-set bounds in Fig. [16] 
we have the following corollary. 

Corollary 2: For a multi-cast, single antenna two-source 
two-destination system, the system DMT d{f) is upper 
bounded by 

d{f) < di3(r), 
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if the system is non-clustered, and by 

d(f) < d22{r), < r < 1, 

if the system is clustered. Here r = ri + r2 is the sum 
multiplexing gain of the system, and this upper bound is 
maximized for f = (r/2,r/2). 

We omit the proof, which is very similar to the upper bound 
calculation in Section [VII- Al 

We observe that cooperative multicast is still limited in 
multiplexing gains. We next study study the cooperative in- 
terference channel, where Di is only required to decode Si 
and D2 to decode 5*2. The cooperative interference channel 
imposes looser decoding requirements on the destinations and 
potentially leads to higher achievable rates. The next corollary 
shows that for the clustered cooperative interference channel it 
is still not possible to achieve multiplexing gains above r = 1. 
Hence, we conclude the multiplexing gain limitation is not 
due to having one source-destination, but is due to the finite 
capacity links within each cluster. 

Corollary 3: A single antenna two-source two-destination, 
clustered cooperative interference channel has the best DMT 
as 

max d{f) = d.22{r), < r < 1. 

ri+r2=r 

Proof: We can show that d{f) < d22{r) using the 
upper bound of Lemma [T] The result in [50] suggests that for 
cooperative interference channel the total multiplexing gain 
can be at most 1. Thus we have d{f) < d22{r), < r < 1. 

A simple achievable scheme assumes that (Si, Di) pair uses 
(S2, D2) pair as relays for half of the transmission period to 
send ^ r/21ogSNR. In the remaining half 5*2 sends 

i?(^s2) ^ r/21ogSNR to D2 utilizing Si and Di as relays. 
Note that equal distribution of rates gives the best network 
diversity, since any other distribution of ri and 7-2 leads to a 
lower diversity for one of the streams. Then, for each case, 
the problem reduces to the one discussed in Section IVII-AI 
We can easily show that this strategy meets the DMT upper 
bound and that using such a time division scheme is DMT 
optimal for the cooperative interference channel. ■ 

For comparison, suppose there were two clustered single 
antenna sources and a single two-antenna destination. This 
system can be a virtual MIMO, achieving the full 2x2 MIMO 
DMT, unlike the single-antenna two-source two-destination 
system described above. In this case, when high diversity gains 
are needed, the sources can cooperate, decode and forward 
each other's signals using time division, and collectively act 
as a two-antenna transmitter similar to the above argument. For 
high multiplexing gains, they simply operate in the multiple 
access mode, i.e. each source sends its own independent 
information stream, and thus can attain all multiplexing gains 
up to 2 [37]. However, if we had a two-antenna source and 
two clustered single antenna destinations, the system would 
be multiplexing gain limited, as CSI is not available at the 
transmitter [62]. All these examples emphasize the difference 
between transmit and receive clusters, in addition to the effect 
of finite capacity links within each cluster 



VIII. Conclusion 

In this work we find the diversity-multiplexing tradeoff 
(DMT) for the following subproblems of a general multiple an- 
tenna network with multiple sources, multiple destinations and 
multiple relays: 1) A single source-destination system, with 
one relay, each node has multiple antennas, 2) The multiple- 
access relay channel with multiple sources, one destination and 
one relay, each node has multiple antennas, 3) A single source- 
destination system with two relays, each node has a single 
antenna, 4) A multiple source-multiple destination system, 
each node has a single antenna. For different configurations 
we consider the effect of half-duplex or full-duplex behavior 
of the relay as well as clustering. 

Firstly, we study a full-duplex multi-antenna relay system 
DMT. We examine the effect of clustering on both the DMT 
upper bounds and achievability results. We compare a single- 
antenna relay system with a multiple-antenna relay system, 
when the source, the destination and the relay have m, n 
and k antennas respectively, and investigate the effects of in- 
creased degrees of freedom on the relaying strategies decode- 
and-forward, and compress-and-forward. We find that multi- 
antenna relay systems have fundamental differences from their 
single-antenna counterparts. Increased degrees of freedom af- 
fects the DMT upper bounds and the performance of different 
relaying strategies leading to some counterintuitive results. 
Although the DF protocol is simple and effective to achieve 
the DMT upper bounds in single antenna relay systems, it can 
be suboptimal for multi-antenna relay systems, even if the 
relay has the same number of antennas as the source. On the 
other hand, the CF strategy is highly robust and achieves the 
DMT upper bounds for all multiplexing gain values for both 
clustered and non-clustered networks. Clustering is essential 
for DF to achieve the DMT upper bound for low multiplexing 
gains, but does not help in the high multiplexing gain region. 
What's more, it has an adverse effect on both the upper 
bound and the DF achievable DMT if the relay has multiple 
antennas due to decreased degrees of freedom in the source- 
relay channel. 

We extend the above full-duplex results obtained for the 
multiple antenna relay channel to the half-duplex relay as 
well. We show that for the multiple-antenna half-duplex relay 
channel the CF protocol achieves the DMT upper bound. 
Although it is hard to find the DMT upper bound explicitly for 
arbitrary to, n and k, we have solutions for special cases. We 
show that the half-duplex DMT bound is tighter than the full- 
duplex bound in general, and CF is DMT optimal for any to, 
n and k. We also argue that the dynamic decode-and-forward 
protocol or any decoding based protocol would be suboptimal 
in the multiple antenna half-duplex relay channel as they are 
suboptimal in the full-duplex case. 

We next investigate the multiple-access relay channel. In 
MARC, CF achieves the upper bound for high multiplexing 
gains, when both users being in outage is the dominant outage 
event. 

Finally, we compare wireless relay and cooperative net- 
works with a physical multi-input multi-output system. We 
show that despite the common belief that the relay or co- 
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operative systems can be virtual MIMO systems, this is 
not possible for all multiplexing gains. Both for relay and 
cooperative systems, even if the nodes are clustered, the 
finite capacity link between nodes in the source cluster and 
the finite capacity link between the nodes in the destination 
cluster are bottlenecks and limit the multiplexing gain of the 
system. Cooperative interference channels are also limited the 
same way. It is straightforward to extend our results for a 
single source-destination pair with multiple relays and for 
cooperative systems with N sources and N destinations with 
each destination decoding all sources. 

Overall, our results indicate the importance of soft infor- 
mation transmission in relay networks, as in CF, and sug- 
gest that protocol design taking into account node locations, 
antenna configurations and transmission/reception constraints 
are essential to harvest diversity and multiplexing gains in 
cooperative systems. 

Appendix I 
Proof of Lemma[T] 

By Proposition [T] the information rates i?*^*^') from node k 
to node I in the network satisfy 

for some p{xi,X2, ...,x„i)- Also, we can easily observe that 
£]y is implied by the event 

Then for any coding scheme with rates i?'^'"'^, we can write 

\kGCi,ieCf 

> P{Ic, < R^^"-'') 

> min P{Ic, < R^'^'='^). 

p{xi,X2,...,Xm) 

The above statement holds for all coding schemes with rates 
/jC-O; thus, it is also true for the one that minimizes the left 
hand side. Then we have 

min P{£n)> min P{Ic, < R'-'^"^^) (38) 

idl coding schemes p(^Xx .X2,--.,3^tji) 

The right hand side is the minimum outage probability for 
cut-set Ci, and by the definition in the lemma 

min P(/c, < i?('^<^'))=SNR-'''^»(''*''''). (39) 

p{xi,X2,...,Xm) 

Using the definition of maximum network diversity, we have 
min P(£w)=SNR-''(*^). (40) 

idl coding schemes 

Substituting ^9} and into (l38]l leads to 

d(f) <dc,(r(C')). 
Since this is true for all the cut-sets, we have 

d{f) < mindc^ir'-'^''^). 

i 

We conclude that the maximum network diversity order d{r) 
is upper bounded by the maximum diversity order of each 1^ ■ 



Appendix II 
Proof of TheoremE] 

In the DF protocol, the source and the relay employ block 
Markov superposition coding [55], [20] and the destination 
does backward decoding [20], [58], [59]. For achievability, 
we constrain the relay to decode the source signal reliably. 
If based on its received SNR the relay cannot decode, then 
it remains silent (or sends a default signal). We assume this 
is known at the destination, which can be communicated at a 
negligible cost. Since fading is constant for all B blocks, this 
has to be communicated only once. 

The relay decodes if the instantaneous mutual information 
satisfies 

p(^^) </(X5;Yfl,|X«). 

If the relay can decode, the mutual information at the des- 
tination is /(XgXij; Yd), otherwise it is /(X5; Yd|X/j). 
We choose X5 and X/j independently as complex Gaussian 
with zero mean with covariance matrices Qs = ImPs/m, and 
Qr = IkPnJk respectively. Then we can write 



/(Xs;Yfl|XH) 
/(Xs;YD|Xfl) 
/(XsXflsY^) 



log Ls.D 
log isi?„D, 



where 



Ls.R — 



m 



and Ls^D and Lsr.d are defined in ( fT3T l and (fl4] i respectively. 
We calculate the probability of outage as 

P (outage at D) 

= P (outage I relay decodes)P(relay decodes) 
+ P (outage I relay cannot decode) 
. P(relay cannot decode) (41) 

= P{LsR.D < SNR'')P(Ps,fl > SNR--) 
+ P{Ls.D < SNBJ')P{Ls,R < SNR'') 



SNR-'^(™+^)"('') + SNR-'*'""('')SNR 



if < r < minj™, n, k} 



if min{TO, n, k} < r < min{m, n} 

for which we used the fact that P{Ls,r > SNR'") = 1 for < 
r < min{m, n, fc}, and for min{rn, n, fc} < r < min{m, n}, 
P{Ls.R > SNR'')=0 and P(Ls,K < SNR'')=1. Hence we 
can write the DMT for DF as in (|24] |. Note that any other 
choice of Qs, Qr. and Q = Cov(X5, X/j) would not improve 
this result. This is because for any Qs, Qr. and Q, due to ([U, 
the mutual information expressions have the upper bounds 



I{Xs;Yr\Xr) 
I{Xs;Yd\Xr) 



< 



< 



log 
log 



Ik 



H 



hIrPs 



sd'H-sdRs 



and 



I(XsXr-Yd) < log KsRM 



where Ksr.d is defined in (|9]l. A DMT calculation using these 
upper bounds would result in the same DMT as in ( |24] |. 
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Appendix III 
Proof of Theorem|4] 

First, the mutual information for cut-set C^i is the same as 
the non-clustered case of (|5]l, and the DMT upper bound for 
this cut-set is (i(m+fc)n(^)- For cut-set Cs we need to find the 



where 



DMT for the channel 
Y 



H 



S,RD 



X + Z, 



where 'H.s,rd is defined in ( fTOl i. is an fc x m matrix with 
all entries equal to VG, and \isD is n x m, with complex 
Gaussian entries hij, i = 1, n, j = 1, 2. The channel input 
X is rn X 1 and has the total transmit power constraint P5. We 
assume Pg = SNR for notational simplicity in this appendix. 
The channel output, Y, and the complex Gaussian noise at the 
output, Z, are (n + k) x 1. 

For m ~ 1, the DMT is easily calculated as rfcs(^) — 
< r < 1, as /cs > log(l + kGPg). 

For m = 2, the instantaneous mutual information for a given 
channel gain matrix Tig is then 



/(X;Y) = log 
Note that 



SNRH^ j^jyils^RD 



' kG kC 
kG kG 



rt 



^SD^SD 



which means having k relay antennas only increases the 
Gaussian channel gain in between the source and the relay 
antennas by a constant factor Therefore, without loss of 
generality we can assume k ~ 1. For A: = 1, 



/(X;Y) 



log (1 + 2GSNR + E^^iE2^i|/i,yf SNR+ 
+ Sr=i|/i^i-/».2pGSNR 
+ Sr=iS;'=»+i|/i»i/^j2 - hahji\^mR^) . 

There is no outage for multiplexing gain r < 1 as 
/(X; Y) > log(l + 2GSNR). For 1 < r < 2, we can lower 
bound the outage probability as 

P(outage) = P(outage|f ) + P(outage|£=)P(f 
> F(outage|f)F(£), 

where 

£ = {mh,i}~^{ha}\<e,Mha}~^{h,2}\<e: 
i = 1, ...,n} 



with e = l/VSNR^. When 8 holds. 



/(X;Y) 



Sr=iS?.il%PSNR 



log (1 + 2GSNR 
+ 2e2nGSNR 

+ 2e2Er=iS;L.+i|/ia - /ijipSNR^) 
Then for a target data rate i?^-^' = r log SNR, we have 

1 



P(outage|£) 



SNR'^-i 

2Er=iE7.,+i|/i.i-/i,iP</(SNR)), 



/(SNR) 



SNR'^ - 1 - 2GSNR - 27iGSNR^ 
SNR'' 



r-l 



Then for SNR > 1, 1/SNR''"^ < 1, therefore, we can further 
lower bound P(outage|£) as 

P(outage|f) > P(l]r=iS,Li|/^.,f 

+ 2Sr=iS7=,+i 1/1.1 - /i.ip < /(SNR)) . 

As /(SNR) = 1, then P(outage|£)^l. On the other hand, as 
the real and imaginary parts of all random variables are i.i.d. 
we have 



and 



P{£)^{Pmh.a}~n{h,2}\<e)y' 



P{\di{h,i}^di{h,2}\<e) 

roc /"X+e 1 

1 _ 2_ 2 

— e ^ ^ dydx 

7T 



— OQ<j x — e 
00 ^ 



-00 -/ — e 

€ /"OO 



e -/ —00 

1 



TT 
1 

— e 

TT 



/27r 



The error function has the Maclaurin series expansion of 



erf (x) 



1 



1 

10' 



1 

42' 



which makes erf(e/V2)=e at high SNR. Then P{£)=e^ 



SNR 



-2n+nr 



and we have 

P(outage)>SNR-2"+"'-. 

On the other hand, 

P(outage) 



< log 



I2 + SNRHI; ^^^H^^^ < r log SNR) 



= SNR 



— 2n-\-nr 



where H^^^,^ 



"-SR "-SD 



, and HsR is an 1 x m matrix 
with i.i.d. complex Gaussian entries. 

Thus, for 1 < r < 2, dcs (r) is equal to the DMT of a 
2 X (n + 1) system, d2{n+i){r) = ri{2 — r), and overall we 
have the d{r) expression stated in Theorem |4] for ra ^ 1,2 
and arbitrary n and k. 

Appendix TV 
Proof of Lemma|4] 

In random state half-duplex relay systems, the system state 
can also be viewed as a channel input. Thus, we need 
to optimize over all joint distributions p{xs,XR,q). Using 
Proposition [U we have 

< min /c • , 
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i — S,D, for some p{xs,xif,q), where R^^-^^ is the in- 
formation rate from S to D. Then for a target data rate 
r(.t) = J. log SNR we have 

min <i?(^) 
all coding schemes ^ 

> min max P{Ici < R'^'^^) 

p(xs,XR,q) i 

> min max min P{ICi < fi^^"*)- 

p(q) i p{p{xs,xn.\q)) 



Then using 



we can write 



d{r) < niayimhidci{r,p{q)). 

p(q) « 



For the multiple antenna, half-duplex relay channel we have 

Ics = I{Xs,Q;Yr,Yd\Xr) 

= I{Xs;Yr,Yd\Xr,Q)+I{Q;Yr,Yd\Xr) 
< I{Xs;Yr,Yd\Xr,Q) + 1 

where the last inequality follows because Q is a binary random 
variable. Similarly, 

Ico = I{Xs,Xr,Q;Yd) 

= I{Xs,Xr;Yd\Q)+I{Q;Yd) 
< I{Xs,Xr;Yd\Q) + 1. 

The above two bounds show that random state protocols can 
at most send one extra bit of information, which does not play 
a role at high SNR. Thus, fixed and random state protocols 
have the same DMT upper bound. 

Appendix V 
Proof of Theorem[6] 

To illustrate that CF achieves the DMT in Theorem |6] we 
follow the CF protocol of Section |IV] In the static half-duplex 
case the relay listens to the source only for t fraction of time 
with < t < 1. The Wyner-Ziv type compression rate is 
such that the compressed signal at the relay can reach the 
destination error- free in the remaining {1 ~ t) fraction of 
time, in which the relay transmits. Then, for a fixed t the 
instantaneous mutual information at the destination is 

= tI{Xs;YRYD\qi) + (1 - t)I{Xs;YD\XR,q2) 

subject to 

tI{YR;YR\YD,qi) < (l-t)/(Xfl,;Yp|g2). (42) 

Note that the above equations incorporate the half-duplex 
constraint into (fTTT i and (fTSl i. The source and relay input 
distributions are independent, Yr is the auxiliary random 
vector which denotes the compressed signal at the relay and 
depends on Yr and X/?. More information on CF can also 
be found in [14], [21], [22] for the half-duplex case for single 
antenna nodes. 

We consider X5 and X^ are i.i.d. complex Gaussian with 
zero mean and covariance matrices ImPg/m, IkPR/k, Yr = 
Yr i + Zr, and Zr is a vector with i.i.d. complex Gaussian 
entries with zero mean and variance Nr that is independent 



from all other random variables. Using the definitions of Ls.d, 
LsR,D, Ls^RD, and L'g (fT3] l. (O, JTSl l and (fT6l l we have 



I{YR;YR\YD,qi) = log 



Ls,RD 

Ls,DN'k 

LsR.D 



Ls., 



Thus using (l42T l we can choose the compression noise variance 
Nr to satisfy 



,r k Ls,RD J- I ^SR,D 

= \l — — — > With U = Ls^D 



U 



and (I42I 1 becomes 



Nr + 1 



S,D- 



(43) 



To prove the DMT of ( |43] | we follow steps similar to (fTsl l- 
Then we have 



P (outage at D) 

pfpC^^) < rlogSNR 



P i L's^RoLto < 2'*SNR' 



P 



{^SR, 



-t) Tt 



< 2** SNR'' 



< P {l'smdL^'^ < 2'SNR 
+ ^(4«,*Xd<2'=SNR^) 

- PitiogL's^Rj, 

+ (1 - i) log Ls,D < r log \/2^SNr) 

+ P{il~t)l0gLsR^D 

+ tlogLs^D < rlog V^SNr) 



(44) 



(45) 



(6) 

= SNR 
= SNR 



SNR 



-min{dc {r,t),d'c (r.t)} 



where (a) is because for any fixed < t < 1, 2'^' > 2*^*. 
For (6) we have used the fact that L'g j^j^ and K'g rjj, Ls^d 
and Ks.D, and Lsr^d and Ksr^d are of the same form 
except for power scaling and hence result in the same DMT. 
As a result if P(outage at D)=SNR-''('''*\ then d{r,t) > 
mm{d'Q^{r,t), d'(2j^{r,t)}. As the achievable DMT cannot be 
larger than the upper bound, we conclude that CF achieves the 
bound in (|32] i for any t. Thus it also achieves the best upper 
bound of ( |33] |. 

If the relay is dynamic, CF can also behave dynamically and 
t will be a function of CSI available at the relay. For dynamic 
CF we can still upper bound the probability of outage at the 
destination with (|45l l, which is equivalent to the DMT upper 
bound for dynamic protocols at high SNR. Hence, dynamic 
CF achieves the dynamic half-duplex DMT upper bound. 
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Appendix VI 
Proof of Theorem[7] 

In this appendix we prove Theorem Q For m = 1, ( [30] l can 
be written as 



where 



with 



I^^ (t) = tlogK's^j,^ + {l-t) log Ks,D, 

I n+k \ 

log 1 + ^ a;,;SNR 



^'s,RD 



Ks,D = log + ^ x.SNRj 

where Xi are independent exponentially distributed random 
variables with parameter 1, that denote the fading power from 
source antenna to receive antenna i at the destination or at the 
relay respectively. 

Let Xi = SNR^"% i — 1, n + k. Then ai are i.i.d. with 
probability density function 

f^^{ai) = log(SNR)SNR-"'exp(-SNR~"'). 

Let A denote the outage event for a target data rate i?^^^ = 
r log SNR. Then probability of outage is 

P{A) 

- P(/^,(t) <rlogSNR) 

fa{oi.)dOL 



(a) 



I (logSNR)"+'^'SNR-^"' exp(-I]SNR-"')da 

J A 

SNR-^"'da 



(6) 



SNR-^^'da 



= SNR-'^ 



where IR("+'^)+ is the set of real (n + fc)-vectors with nonneg- 
ative elements. The outage event A is defined as 

A — {t max{0, 1 — «!, 1 — Q!„+i } 

+ {l-t) max{0, 1 - ai} < r} 

where without loss of generality we assume ai — 

min{ai, ...a„} and a„+i = min{Q!„+i, Q;„+fc}, and G* 
is given as 

G* = inf Sr=i''a*- (46) 

We have (a) because (log SNR)"+'^ does not change the 
diversity gain, exp(— SNR^"*) decays exponentially with 
SNR if a» < 0, exp(-SNR""') is e for a, = and 
exp (-SNR""') approaches 1 for > at high SNR [37], 
(&) follows because at high SNR A converges to A, finally 
(c) is due to Laplace's method [37]. 

As a result dc'^{r,t) — G*. To solve the optimization 
problem of (l46T l. we first solve the subproblems 



inf 



Si 


= {(«1, 




-i)|0 


< 


an 


+1 < 1 < ai} 




= {("1: 


anA 


-i)|0 


< 


ai 


< an+i < 1} 


S3 


= {("1: 


"riH 


-i)|0 


< 


a„ 


+1 < ai < 1} 


Si 


= {(«1: 


"riH 


-i)|0 


< 


ai 


< 1 < a„+i} 


S5 


= {(«1: 






< 


ai. 


: Oin+l}- 



As an example, suppose we want to find si. Thus we have 
the following linear optimization problem 

minimize T."^^ai 
t(l-a„+i)-r < 

< ttn+l < 1 < ai 
minjai, ...«„} — a\ 
min{a„+i, a„+fc} = q:„+i 

This problem has two solutions at 

("i, ■•■,",*, ",*+fc) 

(1,...,1,0,...0) if t<r 

(l,...,l,l-r/t, ...,1-r/i) if t>r ' 



Then for a e Af] M("+'')+ f] Si 
Si = min Stti — 



n if t < r 

n + k{l-r/t) if t>r ' 

Similarly, we find 82,83, 84 and 55. Then G* = mini 8i, which 
concludes the proof. 

Appendix VII 
Proof of Theorem[8] 

When 5*1 and 5*2 do equal time sharing and t = 1/2, we 
use Corollary [T] to conclude that djfj^f^f-, cpi^) ~ ^(1 — r) is 
achievable, where TS denotes time sharing. Next, we discuss 
the case when both sources transmit together 

In the half-duplex MARC, when both sources transmit 
simultaneously and the relay does CF for the signal it receives, 
similar to CF discussed in Sections llVlandlVl the information 
rates satisfy 



< tI{Xs,;YRYD\Xs,,qi) 

+ {l-t)I{Xs,;YD\Xs,Xn,q2) 

< tI{Xs,;YRYD\Xs„qi) 

+ {l-t)I{Xs,;YD\Xs,XR,q2) 

< tI{Xs,Xs,;YRYD\qi) 

+ il-t)IiXs,Xs,;YD\XR,q2) 



for independent Xs^, Xs2, and Xr subject to 

tI{YR-YR\YD,qi) < il-t)I{XR;YD\q2), 



(47) 



where Yr is the auxiliary random variable which denotes the 
quantized signal at the relay and depends on Yr and Xr [44] 
and t is the fraction of time the relay listens. 

To compute these mutual information, we assume and 
Xs2 are independent, complex Gaussian with zero mean, have 
variances Ps^ and Ps^ respectively, and Yr — Yr i + Zr, 
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where Zj^ is a complex Gaussian random variable with zero 
mean and variance Nu and is independent from all other 
random variables. We define 



LSi.D 
Lsi.RD 



L 



LS1S2MD 



^ 1 

^ 1 

^ 1 

^ 1 

^ 1 
+ 

^ 1 



H 



f \hs,D\''Ps, 

f \hs2D\^PS2 

^\hs2D\''Ps^~ 
^\hs^D\''Ps^' 

f Ihs^B^Ps, H 
NRil + \hs,D 
f \hs2B.fPs2 H 

' Ps, 



-\hs2D\''Ps2 

^\hs2D\^PS2 

-\hs,D\'Ps^ 
'Ps.) 
' \hs2D?i 



\hRD\''PR 



S2 



SiS2,RD 



hs2D\^PS2) 



PS2 



Nn 



1 







where 11s^S2,b.d 





" 

1 _ 

hszR 



H 



t 

S1S2MD 



hSiR 

hSiD hs2D 

Since the relay has relevant CSI, using WT\ it can choose 
the compression noise variance Nr to satisfy 



Nr 
Then 



^SiS2,RD 
U 

R{S2) 



with U = LsiS2,D 



S2R,D 



V Ls^S2,D 



< 



< 



< 



t\og 
t\og 
t\og 



^Si,RD _^ 

Nr + 1 

Ls2.RD ^ 
Nr + 1 
LsiS2,RD 



(l-<)logL5,,z3(48) 

(1 - t)logLs2,D 



Nr + 1 
(1 - t) \0gLs,S2,D 



To find a lower bound on the achievable DMT, we use 
the union bound on the probability of outage. For symmetric 
users with individual target data rates iJ^^^i-' = R^'^^^^ = 
r/21ogSNR, and a target sum data rate i?^^' = rlogSNR 
the probability of outage at the destination is 

P (outage at D) 

< < i?(^)/2) + P(i?(^^) < /2) 

+ + < (49) 

One can prove that the first and second terms P{R^^^'> < 
i?('^V2) and P{R^^^^ < R^^'^ /2) are on the order of 
SNR"^^"''/^^ at high SNR, for any t. To see this we write 
(l48T l explicitly as 



< 



t\og 
+ (1 



l + \hs^D\''Ps^ + 



\hs^R.\'Ps, 

Nr + 1 



t) \og{l + \hs,D?Ps:, 



As the relay compresses both sources together, the com- 
pression noise is on the order of SNR and the term 
\hsiR\'^ Psi / R + 1) does not contribute to the overall mutual 
information at high SNR. 



The last term in (|49] | can be analyzed similar to Section HVl 
as this term mimics the 2 antenna source, 1 antenna relay and 
1 antenna destination behavior. For mi — m2 = fc = 1, we 
follow the proof from ( |44] |. 

< P{L's,S2,RDLtslD<'^'^^^'') 

+ ^(4'5*U^ks.,i.<2'SNR'- 

< P (Lks2R.DLtslD < 2'SNR 



PL 



(I-*) 



live 



< 2''SNR' 



^SiS2R,D-^SiS2,D 

= SNR-'^^o^"'*) +SNR-''^i>("'i-*) 
From first line to the second, we used the fact that L[ 



-'SiS2R,D 



with 



SiS2,RD 



H 



S1S2.RD 



Ps, 

PS2 



SiS2,RD 



'^SiS2.RD + I2 



> 



as a 2 X 2 multiple antenna system has higher capacity than a 
3x1 system. 

Using dQ^{r,t) from Theorem |7] with m = mi + m2 — 2, 
k = 1, n — 1, to maximize minjii^^ (r, t), d'^^ (r, 1 — t)} over 
t, we need to choose < t < | and thus 



,SIM 

^MARC.CF 



(r) > min |l — — , 3(1 — r)| 



where SIAI denotes simultaneous transmisssion. 

To find an upper bound on the achievable DMT we write 

P (outage at D) 

> max{P(p(^i) < /2), P(P(^^) < PC^) /2)}, 



so 



'^MARCCFV) ^ ^ — 



Combining this with the upper bound in (l34l i. and with 



ITS 

A.J I 



MARCCF 



(r), we have Theorem [8] 



Appendix VIII 
Proof of Upper Bound in Theorems [9] and [Tol 

To provide upper bounds, we will use the cut-set bounds 
as argued in Lemma [T] The cut-sets of interest are shown in 
Fig. [T3] and denoted as Cg, Csr^ and Cd- We will see that 
these will be adequate to provide a tight bound. 

In order to calculate the diversity orders dci{r) for each 
cut-set, we write down the instantaneous mutual information 
expressions given the fading levels as 

Ics = H^s ; ^R.! Yr^ Yd I Xr^ Xr^ ) 
IcsRi = I{XsXr^;Yr.^Yd\Xr^) 
Icd = I{XsXr^Xr.^\Yd)- 

To maximize this upper bound we need to choose Xs, Xr^ 
and complex Gaussian with zero mean and variances P5, 
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Pfjj and respectively, where P5, Pr^ and Pr^ denote the 
average power constraints each node has [57]. Then 



ICs < I'cs 



log (1 + \asR, pPs + \hsR, ?Ps + \hsD\''Ps) 

(50) 



< I'c 



log 



I2 + lisRi.R2D'H-SB,uR2DiPs + PRi) 



(51) 



= log {1 + {\hsD\^ + \hR,D\^ + \aR,D\^) 

■ {Ps + Pr, + Pr^)) (52) 



where 



H 



SR1.R2D 



hsR.2 

hsD 



hRiR2 

hR^D 



(53) 



and we used ^ to upper bound Icsr^ with /^^^ in (ISTT i and 
Ic^ with /^^ in (|52li. 

For a tax-get data rate R'^'^^ = rlogSNR, P(/cs < 
< i?(^))=SNR-''"(''\ whereas P(/cs^^ < 
p(^))^SNR"''^^(''). Then using Lemma[I] the best achievable 
diversity d{r) of a non-clustered system is upper bounded by 
d{r) < niin{di3(r),d22(r)} = di^ir). 

When the system is clustered, /^^ and /^^ are larger 
than the Gaussian channel capacities log (1 + Gsi^iPs) 
and log (l + GR2DPR2) respectively. Then P(/^ 



Cs 



< 



R(T)-^=P{I^^ < i?('^))=SNR"°°, if r < 1. In other words, 
it is possible to operate at the positive rate of i?'^) reliably 
without any outage and as SNR increases, the data rate of 
this bound can increase as log SNR without any penalty in 
reliability. However, this is not the case for any r > 1 as 
P(/^^ =< P('^))=P(/^^ < P('^))=l. Combining these 
results with the upper bound due to Icsr-^ , we have 



d{r) < 



d22{r) 




if 
if 



r < 1 
r > 1 



Appendix IX 
Proof of Achievability in Theorem|9] 

We assume the source. Pi and R2 perform block Markov 
superposition coding. After each block Pi and P2 attempt to 
decode the source. The destination does backward decoding 
similar to the case in Section |IV] 

Using the block Markov coding structure both Pi and P2 
can remove each other's signal from its own received signal 
in (|35] | and (|36] | before trying to decode any information. 

We choose Xr-^ and Xr,^ independent complex Gaussian 
with zero mean and variances Pr^ and Pr^ respectively. 
We also choose Xs independently with complex Gaussian 
distribution CJ\f{0, Pg). Then the probability of outage for this 
system, when the target data rate P'-^' — r log SNR, is equal 



to 

P (outage at D) 

= P (outage and both relays decode) 

+ P(outage and Pi decodes, P2 cannot decode) 
+ P(outage and Pi cannot decode, and P2 decodes) 
+ P(outage, and none of the relays decode). (54) 

P (outage at D) 

= P{LsR,R2,D < SNR^ Ls,R, > SNR^ Lsm2 > SNR"-) 
+ P{LsR,M < SNR^ Ps,i?i > SNR^Ps,i^. < SNR'') 
+ P{LsR2:D < SNR^ Ps,i?i < SNR\LsM2 > SNR'') 
+ P{Ls,D < SNR^ Ls.iix < SNR^Ps.i^. < SNR"^), 

where 

Ls,R, = log {l + \asRA^Ps) (55) 
Ls,R, ^ log{l + \hsRj^Ps) 

LSR^R2,D = l0g{l + \hsD?Ps + \hR,D?PR, 

+ \aR.,D?PR.2) (56) 

LSR,.D - l0g{l + \hsD\''Ps + \hR,D\''PR^) {51) 
LSR.2.D = log{l + \hsD\''Ps + \hR2D\''PR2) 
Ls.D - log{l + \hsD\''Ps), 

and asRi — hs.Ri, CIR2D — ^r^d as the system is non- 
clustered. 

Using the fact that P{Ls,r, > SNR'')=1 and P{Ls.r2 > 
SNR'")=1, this outage probability becomes 

P (outage at D) 

= SNR^'^^i^'') + 2SNR-'^^i('')SNR-''"('') 

+ SNR-'^" SNR-'^" ('■^SNR-''" 
= SNR-''^i('') 

at high SNR, which is equivalent to the outage behavior of 
a 1 X 3 system (or 3x1) system. Hence, in a non-clustered 
system if both relays do DF, the DMT in Theorem |9] can be 
achieved. 

Appendix X 
Proof of Achievability in Theorem[TO] 

To prove that the DMT of Theorem[TO]is achievable, we use 
the mixed strategy suggested in [20], in which Pi does DF and 
then the source node and the first relay together perform block 
Markov superposition encoding. Similar to the non-clustered 
case in Appendix |IX] we require Pi to decode the source 
message reliably, and to transmit only if this is the case. We 
assume that P2 and the destination know if Pi transmits or 
not. The second relay P2 does CF. 

To prove the DMT we calculate the probability of outage 

as 

P(outage at D) — P(outage|Pi decodes)P(Pi decodes) 
+ P(outage|Pi cannot decode) 
. P(Pi cannot decode). 
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As i?i is clustered with the source, source to Ri com- 
munication is reHable for all multiplexing gains up to 1; 
i.e. P{Ri decodes) and P{Ri cannot decode) can be made 
arbitrarily close to 1 and respectively. Therefore, we only 
need to show that P(outage|i?i decodes) decays at least as 
fast as d22{r) with increasing SNR. 

When i?i decodes the source message reliably, if the target 
data rate i?^"^^ satisfies 



subject to 



I{Yr, ; Yr, \Xr,Yd) < I{Xr, ; Yd) 



(58) 



(59) 



then the system is not in outage. 

We choose Xs, Xr^ and Xr^ independent complex Gaus- 
sian with variances Ps, Pr^ and Pr,^ respectively and Yr,^ = 
Yr2 + Zr^, where Zr^ is an independent complex Gaussian 
random variable with zero mean, variance Nr^ and indepen- 
dent from all other random variables. 

We define 



L 



SRi,R2D 



H 



SRi,R2D 



SRi,R2D 



H 



Ps 
Pr, 

Nr,+1 o" 
1 

' Ps 
Pr, 



SRi,R2D 



"-SRi,R2D 



'-SRi,R2D 



where flsRi,R2D is given in i53i . Using the definitions of 
Ls,Ri, LsRi,D, and Lsr.r^^d from (jSTJi and §6^, 

with as.Ri ~ G and ur^.d ~ G, the instantaneous mutual 
information expressions conditioned on the fading levels for 
the mixed strategy become 

I{Xs;Yr,\Xr,) = log Ls,R, 

I{XsXr,;Yr,Yd\Xr,) - log^^^^ifi^. 

The mutual information in the compression rate constraint of 



I{Yr,;YrJXr,Yd) = log 
IiXR,;YD) = log 



^SRi,R2D 
LsRi,dNr2 
LsRiR2,D 



LsRi.D 

Then the compression noise power has to be chosen to satisfy 

LsRi,R2D 



N 



R2 



> 



L 



(60) 



SRiR2,D 



Note that both sides of the above inequality are functions 
of Nr^. Using the CSI, the relay will always ensure ( l60l l is 
satisfied. 

After substituting the value of the compression noise in 
(|58] | we need to calculate P(outage|i?i decodes). Given Ri 
decodes, this problem becomes similar to Problem 1, and 
we can find that when P(outage|i?i decodes)^SNR"'^'^'"'', 
d{r) = d22{r). Finally, as P(outage at D)=SNR"''('^', we say 
the mixed strategy achieves the DMT bound. 
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